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Abstract 

sh Gaudin model with jordanian twist is studied. This system can be obtained as 
the semiclassical limit of the XXX spin chain deformed by the jordanian twist. The 
appropriate creation operators that yield the Bethe states of the Gaudin model and 
consequently its spectrum are defined. Their commutation relations with the generators 
of the corresponding loop algebra as well as with the generating function of integrals 
of motion are given. The inner products and norms of Bethe states and the relation to 
the solutions of the Knizhnik-Zamolodchikov equations are discussed. 



1 Introduction 



The Quantum Inverse Scattering Method (QISM) was largely created by L. D. Faddeev and 
his school at St. Petersburg as a quantum counterpart of the Classical Inverse Scattering 
Method 0-0- Classifying quantum solvable systems with respect to the underlying dy- 
namical symmetry algebras, one could say that the Gaudin models E] can be seen as 
the simplest ones being based on loop algebras. Their hamiltonians are related to classical 
r-matrices, 

AT 

H(a) = $>«&(*» " (1-1) 
byia 

The condition of their commutativity [H^ a \ H^] = is nothing clsc but the classical Yang- 
Baxter equation 

[r a b(z a - Zb), r ac (z a - z c ) + r bc (z b - z c )} + [r ac (z a - z c ), r bc (z b - z c )} = (1.2) 

where r is antisymmetric and belongs to the tensor product q (g) g of a Lie algebra g, or its 
representations and the indices fbc the corresponding factors in the iV-fold tensor product 
of this algebra. The Gaudin models based on the classical r-matrices of simple Lie algebras 
attracted a lot of attention [7)-|ll|. Their spectrum and corresponding cigcnfunctions were 
obtained using different methods such as coordinate and algebraic Bethe ansatz, separated 
variables, etc. The correlation functions of the SI2 Gaudin system were calculated by means 
of the Gauss factorization [H]. A connection between the Bethe vectors of the Gaudin 
models for simple Lie algebras to the solutions of the Knizhnik-Zamolodchikov equation 
was established in 0-^3. An explanation of this connection based on Wakimoto modules 
at critical level of the underlying affine algebra was given in |l()j . 

The algebraic Bethe Ansatz for the Gaudin model based on the invariant classical 
r-matrix deformed by the constant jordanian r-matrix was postulated in (T2J - Following the 
ideas used in the case of the osp(l|2) trigonometric Gaudin model Kulish noticed that 
the similarity transformation by exp(aX + ) ®exp(aJT + ) on the 5(2 trigonometrical classical 
r-matrix 

e ~ A r (+) + e r f-1 
sinh(A) sinh(A) 

. e , ... ( \h ® h + 2X+ (g) X~) 
smh(A) \2 / 

-r^—- (±h®h + 2X-®X+) , 
smh(A) \2 / 

setting A — > eX, a —y and after the scaling limit 

lim er tng (e\) = ^(h(g>h + 2{X+ ® X- + X~ <g> X+)) + f (h ® X+ - X+ ® h), 

yields the s^-invariant classical r-matrix deformed by the jordanian r-matrix. Moreover 
the highest weight vector of the corresponding Gaudin model is preserved. Based on these 
arguments Kulish postulated the Bethe vectors, the spectrum and the Bethe equations of 
the system. 

Alternatively, the jordanian twist |14H15lH^| 

T J = e ma = exp(/i (g) ~ ln (1 + 20X + )) (1.3) 



1"trigW — 



+ 
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can be applied to the s^-invariant spin system based on the Yang quantum R-matrix 



R(X; V ) = I+^P (1.4) 

where V is the permutation matrix in the tensor product C 2 ® C 2 , to obtain the twisted s [2 
spin system El ■ The semiclassical limit r\ — > of this system yields the same Gaudin 
model discussed by Kulish [T5|. The twisted XXX spin chain related to the quantum R- 
matrix 

R(\ ;m 9)=R J (9) + ^T, (1.5) 

where R J (9) = ^^.(J 7 ^) - 1 is the sfe jordanian quantum R-matrix studied in |181 1191 120) . 
whose homogeneous case was analyzed in |14j , will be discussed elsewhere |17| , here will only 
be presented the aspects relevant to the study of the corresponding Gaudin model. 

The £-operator of the quantum spin system on a one-dimensional lattice with N sites 
coincides with i?-matrix acting on a tensor product Vq (8) V a of auxiliary space Vq — C 2 and 
the space of states at site a = 1,2, ... N 

w -*>-(7 T>i^(| + %)• m 

where z a is a parameter of inhomogeneity (site dependence). Corresponding monodromy 
matrix T is an ordered product of the £-operators 



N 



%{\{z a }?) = Z 0N {\- ZN )...Z QX {\-z 1 ) = J[Z (ia {\-z a )= ( ^n) d{\) j" 

a—l ^ ' 

The commutation relations of the X-matrix entries follow form the FRT-relation 

R 12 (X - w m 0)X(A)I(Ai) = T( At )T(A)i? 12 (A - w ry; 0). (1.8) 

Multiplying (|1.8(1 by ii^f and taking the trace over C 2 (g) C 2 , one gets commutativity of the 
transfer matrix 

t(A) = tr X(A) = A(X) + D(X) (1.9) 

for different values of the spectral parameter t(A)t(/x) = t(/i)t(A). 

It is of interest to choose different spins l a at different sites of the lattice, hence the 
following space of states 

n= ® , 

a—l 

with the highest spin vector fi + = <S>a = iU> a . It is straightforward to show that 

C(A)ft+=0, (1.10) 

and 

A - z a 



A(\)n + = a (m + = [] ( X Za + lar)/2 ) n + , (1.11) 



a=l 



£>(A)0+ = d(A)Q+ = J] ( A a - / 07?/2 N ) ° + " (L12) 



0=1 
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As the first step in the application of the algebraic Bethe Ansatz to the twisted XXX spin 
chain one can confirm that the highcst spin vector f2 + is an eigenvector of the transfer matrix 
(Ol) 

t(A)fi+ = A (A)fi+, (1.13) 

with Ao(A) = a(A) + d(X). The next step is to show that = B(fi)ft,+ is also an 

eigenvector of the transfer matrix 



t(A)*i(/x) = t(X)B(n)n + = A!(A,/i)*iOi) + 



A — /i 

- 9{a{n) - d(n))(a{X) - d(A))0 + 
whcre the eigenvalue Ai(A,/i) and the Bethe equations are given by 



(a( M )-d(^))* 1 (A) 



A 1 (A,/ 1 )=A (A) 
a(/x) — d(p,) = 0. 



A — /i 



(a(A) - d(A)), 



(1.14) 

(1.15) 
(1.16) 



The following steps of the algebraic Bethe Ansatz are analogous having the same eigenvalues 
and Bethe equations as in the invariant case |14j . the details will be presented elsewhere |17| . 

The corresponding Gaudin model can be obtained from the spin system via semiclassical 
limit |71l21| by setting = — 2£ and using the expansion in powers of r\ of the monodromy 
matrix 



T(A)=I+|i(A) + 0(r ? 2 ), 



where 



L(A) = 



f h(X) 2X-(X) 
^2A+(A) -h(X) 



(1.17) 
(1.18) 



and 



N 



h(X)=J2 



X - Z n 



N 



tx+),x-(X) = Yl 



X — Zn 



N 



Substituting the expansion of the monodromy matrbc (| 1 . 1 Tfl and 



R(X;r);6) 



= I+l ri (X) + 0(r, 2 ), 



(1.19) 
(1.20) 



with 



(1.21) 



~7 K e 

o f -i -e ' 
\o o o i y 

into the FRT relations 1)1. 8|l . the first nontrivial term, the coeficient of j? 2 , is the so-called 
Sklyanin bracket 



£(A),M 



r e (A-M),i(A) + L(/i) 



(1.22) 



Moreover the central element A(A) [2], of the algebra (|1.8|l . admits the expansion in ry such 
that the second order term of t(A) — A(A) is the generating function of the integrals of motion 
of the corresponding Gaudin model [Jj 



t(X) = -Tr L 2 (A). 



(1.23) 
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The first step in the application of the algebraic Bethe Ansatz, to the Gaudin models, 
is to define appropriate creation operators that yield the Bethe states and consequently the 
spectrum of the generating function i(A). The creation operators used in the s[2-invariant 
Gaudin model coincide with one of the L-matrix entry Ej. However, in the present 
case these operators are not homogeneous polynomials of the generator X~(A) and can 
be defined by the recursive relation as proposed by Kulish |12) . It is convenient to define 

Bj^(/j,i, . . . , /ijw), a more general set of operators, in order to simplify the calculation of 
the commutators between the creation operators and the generating function of the inte- 
grals of motion t(X). These operators are symmetric functions of their arguments and they 
satisfy certain recursive relations. Their commutation relations with the generators of the 
loop algebra are straightforward to calculate and they are essential in the main step of the 
algebraic Bethe Ansatz. The creation operators of the 5(2 Gaudin model with jordanian 
twist are the particular case £?m(mi, ■ ■ ■ , Mm) = -^m (MIj • • • >MAf)- Thus the commutation 
relations between the creation operators and i(A) are easily calculated. Therefore the corre- 
sponding Bethe vectors can be defined by the action of the creation operators on the highest 
spin vector and thus the spectrum of the system determined. In this way the algebraic 
Bethe Ansatz is fully implemented confirming the result of Kulish J2| that the spectrum of 
the system coincides with the one of the s^-invariant model and consequently the twisted 
Gaudin hamiltoneans have the same spectrum as in the s^-invariant case, although the 
Bethe vectors of the two systems are different. However the Bethe vectors, in this case, are 
not eigenstates of the global Cartan generator h g i. 

Besides the problem of determining the spectrum of the system, via algebraic Bethe 
Ansatz, some properties of the creation operators are fundamental in calculating the inner 
products and the norms of the Bethe states. In particular, it turns out that the relation 
between the creation operators of the system and the ones in the untwisted case is essencial 
in determining the inner products of the Bethe vectors. In addition, it is necessary to 
consider the dual creation operators B^(/ji, . . . , hm) obtained by using the dual Gaudin 
model based on the following classical r-matrix 

r* (A) = T (g) h + 2(X+ <g> X- + X- ® X+)) + £(h ® X~ - X~ ® h). (1.24) 
A 

As opposed to the sl2-invariant case, the Bethe vectors ^m(pi, ■ ■ ■ , Mm) are not orthogonal 
for different M's. It should also be mentioned that, due to the jordanian twist, the dual 
generating function of integrals of motion t* (A) is not equal to t(A), the generating function 
of integrals of motion of the original model, contrary to the sl2-invariant case, thus these 
operators are not hermitian. 

A connection between the Bethe vectors, when the Bethe equations are not imposed 
on their parameters, of the twisted Gaudin model to the solutions of the corresponding 
Knizhnik -Zamolodchikov equation, similarly to the s^-invariant model, is based on some 
analytical properties of the creation operators Bm(mi> ■ ■ ■ >Mm)- 

The article is organized as follows. In section 2 we discuss the Gaudin model based 
on the SI2 invariant r-matrix with jordanian twist emphasizing the creation operators Bm- 
Using the previously established properties of the creation operators -Bm, in section 3 the 
spectrum and the Bethe vectors of the system are given. The dual creation operators are 
used to obtain the expressions for inner products and norms of Bethe states. In conclusions 
a relation between the Bethe vectors and the solutions to the Knizhnik-Zamolodchikov 
equation are discussed. In the appendbc the proofs of the lemmas are given. 
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2 5 12 Twisted Gaudin model 

The sl2-invariant classical r-matrix with jordanian twist |22U18ll*H?| is the following element 
of sl 2 <S>Sl2 



1 



r d x ) = "f + £ r =-?(h®h + 2{X+ ® X- + X~ <g> *+)) + £(h (g) X" 



/*). (2.1) 



The matrix form of r^(A) in the fundamental representation of SI2 follows from (|2.1|l by 
replacing the appropriate matrices for the generators of SI2 is given explicitly by 



? 
o f 

V o o 



A 



o 



-e 



(2.2) 



here A e C is the spectral and £ e C is the twisting parameter. Dcfinition of the Gaudin 
model requires not only the classical r-matrix but also the L-operator 



f h(X) 2X-(X) 
MAj - ( v 2X+(A) -h(X) 



(2.3) 



whose entries (h(X), X ± (X))) are generators of the loop algebra £^(5(2) defined by the 
Sklyanin linear bracket 



r 12 (A-/z),L(A)+L(/z) 



(2.4) 



The corresponding commutation relations between the generators are explicitly given by 

[h(X),h(n)]=2£(X+(X)-X+(n)) 

[x-(x),x-(jm)] =-t(x-(\)-x-( Ji j), 

h(X) - h(jt) 



[x+{\),x-(ri\ = - 

[X+(X),X+(n)]=Q 
[h(\),X-(jij\ 

[h(\),X + 0*j\ 



= 2 



= -2 



A — /i 

x-(x)-x-(v) | 

A — (i 
X+(X)-X+( f i) 



cx+(x\ 



(2.5) 



X — /i 



The usual 5(2 loop algebra is recovered be setting £ = 0. 

In order to define a dynamical system besides the algebra of observables /^(sfo) an 
hamiltonian should be specified. Due to the Sklyanin linear bracket (12.4(1 the clcments 



t(X) = - Tr L 2 {X) = h 2 (X) + 2(X+(X)X-{X)+X-{X)X + (X)) 



h 2 {X) - 2ti(X) + 2 {2X-(X) + f) X+(X) 



(2.6) 



generate an abelian subalgebra 



t(X)t(fi)=t(fi)t(X), 



(2.7) 
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One way to show l|2.7[) is to notice that the commutation relation between t(X) and 
can be written in the form 

[t(\),L(ji)] = [M(\,iji),L(jji)], (2.8) 
using (2.4-6), it is straightforward to calculate M(X,p) 

M(A,/i) = -Tr (ri 2 (A - /i)i(A)) - ±Tr (rf 2 (A - M )) 

_/-2g + 2^(A) -4*ffl-2fli(A)\ 3 

~ V -4^f 2^-2£X+(A)yl m) 2 ^ 

Thus the commutators between the generating function i(A) and the generators of the loop 
algebra follow from (2.8-9) 

[mitw]=t »»_ t{It(rw , (2 ,0) 

A — /i 
X — fi 

rn , v , , n X-(X)h(u) - X-(u)h(X) . 

[t(A), - 4 v ; + 2^(A)/i(/x) 12) 

+ ^X-(p)X + {X) + 2ex+{ l i). 

Preserving some generality, the representation space H of the Sklyanin algebra (2.4-5) 
can be a highest spin p(X) representation with the highest spin vector il + 

x+(x)n+ = o, h{x)n+ = P (x)n+. (2.13) 

The spectrum and the eigenstates of t(X) can be studied in this general setting. However to 
have a physical interpretation a local realization 

H = V 1 ®---®V N . (2.14) 

as tensor product of s^-modules is nccdcd. Thcn the L-operator is given by 



TV 



were Y a = (h a ,X^) E End(V" a ) are 5I2 generators in representation V a , associated with each 
site a. For convenience, the generators Y(X) = (h(X), X ± (X)) are written down explicitly 

MA) -£ ( + «) , x- m = £ - f „.) , Z +W - £ j5L. 

a—l v 7 a—l v 7 a— 1 

(2.16) 

In this realization it is useful to consider the expressions of the generators of the global SI2 
Lie algebra in terms of the local ones 

N 



Y gl =Y,y ai (2.17) 
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whcrc Y = (h,X ± ). Also, the following notation is useful 

F(A)o=y(A)| £=0 , (2.18) 

then h(X) = h(A) + £X+. 

A representation of the model with the above Gaudin realization is obtained by con- 
sidering at each site a an irreducible representation Va of s [2 with highest weight £ a 
corresponding to a singular vector ui a £ Va such that Xjui a = and h a uj a = £ a ^a- Thus 
the space of states is 

H = V^ l) ®-..® vff N \ (2.19) 
with the highest spin vector (|2.13|l 

il + = u>i ® ■ ■ ■ ® u>n (2.20) 

and the corresponding highest spin 

N f 

PW = E^7- ( 2 - 21 ) 



The Gaudin hamiltoneans 

N N 



= ]T r ab (z a "**) = £( + C (h a X+ - X+h b ) ) , (2.22) 



where c® (a, b) = h a (§) + 2(X+ (g) X h + X a (g) can be obtained as the residues of the 
generating function i(A) at the points A = z a ,a = l,...,iV using the expansion 



a— 1 ^ ^ a ' 



a,6=l 



here C2(a) = h\ + 2h a + 4X~X+ is the SI2 Casimir at site a. As opposed to the s^-invariant 
case, the generating function 12.23JI commutes only with the generator X^. 

The first step in the algebraic Bethe Ansatz is to define appropriate creation operators 
that yield the Bethe states and consequently the spectrum of the generating function t(X). 
The creation operators used in the s^-invariant Gaudin model coincide with one of the 
L-matrix entry |7j . However, in the present case these operators are not homogeneous 
polynomials of the generator X~ (A) . It is convenient to define a more general set of operators 
in order to simplify the prescntation. 

Definition 2.1. Given two integers M and k > consider the operators 
B&Vi. . . . , M m) = (x-(m) + (X-(jt2) + (k + ■ ■ ■ (X-(ji M ) + (M + k- 1)C) 

M+k-1 

= f] (X-0*n_fc+i)+n£), (2.24) 



n— k 



with B { Q k) = 1 and B { $ = Q for M < 0. 

The following lemma describes some properties of the Bw ■ ■ ■ , /ijw) operators which 
will be used later 011. 
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Lemma 2.1. Some useful properties of B^i^i, ■ ■ ■ , Mm) operators are 

i. The operators B^Qx±, . . . , Mm) are symmetric functions of their arguments. 

n. sj£V» . . . , mm) = sSLiOii, ■ ■ . ,mm-i) (fifW) + (m + fc - i)e) . 

iii. B^{fi) = B^~ 1 ^(jjl) + £,J2iLi ^m-i(m^)> ^ ere A* * s a °/-^ = ImI complex scalars 

with a particular ordering fi — {fJ,i, . . . , MAf} assumed for convenience and the notation 

f' ' /< {/'-: 

for any distinct i\, . . . E {1, . . . , M}. 

iv. <(,.) = BfViKSV 13 ). 

Implementation of the algebraic Bethe Ansatz requires the commutation relations 
between the generators of the loop algebra C^sls) and the Bj[j(hi, . . . , mm) operators. 
To this end we first introduce the notation for the Bethe operators 

M (\ ;tl )=h(\) + Y / -?—, (2.25) 
ven r 

where fi = {mi, . . . , MM-i} an( i A e C \ /z. In the particular case M — 1, /3i(A; 0) = h(X) 
and will be denoted by /?i(A). The required commutators are given in the following lemma. 

Lemma 2.2. The commutation relation between the generators h(X), X ± (A) and the 
M ( Ml 1 ■ • ■ t t 1 m) operators are given by 

A/ 

+ «E fl ri(M (<1 )^(w; (2.26) 

i=l 

i=1 A - W 

M R (*+i)crii i i M 
-2 E a- \ +^ B M + -h^)X + {^ (2.27) 



M 



X-(X)B^ (n) = flW ( M) ( X - (A) + MC) - C E B m ( A u M W ) 

»=i 

M 

= B<*> +1 (A U M ) - e E B m( A u M W )- ( 2 - 28 ) 



The £?-operators that define the Bethe states of the system were proposed by 
Kulish These operators are the particular case k = of 12.24fl and will be denoted by 
Sm(mij ■ ■ ■ 1 Mm) = -B^/ (jtii, . . . , Mm)- A recursive relation of the B-operators follows from 
ii. of the lemma |2~D 

Bw(ft) = B M -i(M (M) ) (^"(Mm) + (M - 1)£) . (2.29) 
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It may be useful to write down explicitly first few B-operators 

B = l, B 1 (ji.)=X-(n), S 2 ( Ml , M2 ) =X-(ji 1 )X-(it 2 )+tX-0*i) 



B 3 fa 1 ,to,ii 3 )=X-fai)X-(jt 2 )X-fa 3 ) + 2ZX-(j il )X-(ii 2 ) 
+ ZX-(p 1 )X-(^) + 2eX-(n 1 ), 

As a particular case of lemma l2~2*l the commutation relation between loop algebra generators 
h(X), X ± (X) and the B-operators are given by settting k = iri (2.26-28) 

h(X)B M („) = B u (p)h{\) + J2 f 2 BM(XU^r)-BM + {B W ^W)^. „Wj 

(2.30) 

X+(A)5 M ( M) = BmWI + (A) - f B« >«) ^ M(A; ^:ff ( ^ ;/AW) - e* + (A)) 

- 2 E n Jx 2 - 31 

^ (A-/i l )(A-^) 



i=l 

The crucial step in the algebraic Bethe Ansatz is to determine the commutation relations 
between the creation operators, in this case the B-operators, and the generating function 
i(A). The main lemma, based on the results established previously in this section, will given 
a complcte cxpression for the required commutator. 

Lemma 2.3. The generating function t(X) has the following commutation relation with the 
B-operators: 

( M AhtW M 

t{\)B M {n) = B M {n) t(X) - ]T 7^ + E 7Y— 



i<j 



D /Ali i)^ 

4E T_ ^ M(^;/x (i) )+2^<-i^ W )MA)/3M^;/xW) 

i— 1 2—1 

4 ^ E a^I /MWM ( >) (2.33) 

M 



+ e E 41 2 (m (ij) )/3m-i(^;m (ij) )/3m(^;a* w ) 

M 

+ AM^B M (p)X+ (A) + 2£ 2 ]T ) A+ 

i=i 
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Prooj. The case M = 1 can be obtained directly from the M-matrix (|2.9H and is given 
by H2.12f> . For M > 1 the commutator between the operator t(X) and the corresponding 
B-operator is to be calculated directly using expression l|2.6[) 

[t(X), B M {n)} = [h(X), [h(X), B M (n)]] + 2 [h(X), B M (n)} h(X) - 2^ [h(\),B M (n)] 
+ 2(2X-(X)+^) [X+{X),B M (n)] +4[X-(X),B M (fi)] X+(X). 

The terms in the previous expression only involve the commutators between the generators 
of the loop algebra and the B-operators. Each term in the above equation will be calculated 
separately. Using the commutators (|2.2(i|) and (2.30-32) the first term is given by 

M 

[MA), [h(X),B M (fJ-)]] = B M (fj,) 7T wr r 



4E 



' B M ({X} U m«) - B m {l^) £B M ({X} U ^ 
(X - fii) 2 X - Hi 



^ X-(X)B$_ 1 ({X}U^) ^ B M ({X}U^) ^ 2 

+ 8 Z. n- ll A(\-, l A 4 L — nrTZ — z. 



i<3 (A-Mi)(A-Mj) ^ A- Ml £^A-Mj 

<t f 4i(Wu^) , ot f B d) (iiH) J m (x-,^)-Pm(^) 

+ (A- /*)(>-/*,) +2 ^^-r^ ) 



A/ 



i=l 



The commutator between the generator /i(A) and Bm(p) is also used to determine the next 
two terms 



2 [MA), BmW] MA) = B M f) +4^ BM f U/J '\ (A) 

^ A -Mi ^ A- M; 

a/ a/ 

+ 2^ < 3 _ 1 (^ )h(X)$ M ( Mi ; M W ) " 4? 2 E B m- i (M W ) (* + ( A) - X+ (m, ) ) , 

i=l i=l 

and after diferentiating l|2.30|l 

2 d f , mB M1 f ^({AIUm'VBmM il^({A}U^) 
-2-[MA),S M (M)]=4^ — — — . 
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The next term follows directly from l|2.31[) 

2(2X-(A) + {) [X+(A), Bm(m)] = -4^m({A} U „W) ^ ^ = ^ ^ 



i=l 

M 



- X-(X)B^_ 1 ({X}U^) " (1) ( , th /3 M (A; M W)-/3 M (^;M W ) 

^ ( — Mi) (A — Mj) 

M 

i=l 

The last term is obtained from (|2.32[1 

M 

4[X-(X),B M (n)]X+(X)=4 M C B M (n)X+ (A) - 4 £ £ B M ({A} U M « )X+ (A). 



When all terms in [i(A), -Bm(m)] are P u t together it is straightforward to obtain the formula 
(l2~m □ 

Besides the problem of determining the spectrum of the system, via algebraic Bethe 
Ansatz, some properties of the B-operators are fundamental in calculating the inner products 
and the norms of the Bethe states. To this end the relation between the B-operators and 
the untwisted ones is established using the Gaudin realization Ij2.16|l : 

M-l M 

B M {fj)= E B M ^ n ,...,^ M _ k ) p{ M - k) +e l PM (2.34) 

fe=o h<-<j M _ h 

where 

B M -k(Vji,- ■ ■ > Mjm- Jo = B M _ k (iij x ,. . . ,/J.j M _ k )\s=o 

and also 

# = (-^+i)---(-^+i + i-l) with p i= pf\ (2.35) 

are operators defined for any integers i and j. 

Moreover, some analytical properties of the B-operators are important when solving the 
corresponding Knizhnik-Zamolodchikov equation. 

Lemma 2.4. The B-operators satisfy the following identity 

r) M 8 

— B M (/x) = -^Q^(xMB$-M i) )) (2-36) 

where X~(im) = — ^ - §/i a for site a = 1, . . . , N. 

Using the relevant properties of the B-operators established in this section, the funda- 
mental description of the 5I2 Gaudin model with jordanian twist can be obtained. 
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3 Spectrum and Bethe vectors of the twisted 5^ Gaudin 
model 



In this section the spectrum and the Bethe vectors of the twisted SI2 Gaudin model as 
well as their inner products and norms will be determined by applying the algebraic Bethe 
Ansatz. The first step is to defmc the Bethe vectors 'P a/ (/j) = Bm(h)Q+ by the action of 
the f?-operators on the highest spin vector f2 + . Then the key observation is that highest spin 
vector J7 + is an eigenvector of the generating function of integrals of motion t (A) . Finally the 
spectrum of the system is obtained as a consequence of the commutation relations between 
t(X) and B m (p) lemma l2~3fl H2.33f) . The unwanted terms coming from the commutator are 
annihilated by the Bethe equations on the parameters fi = {fJ,i, ■ ■ ■ ,A*m} as well as by the 
condition X + (\)fl + = 0- Hence the algebraic Bethe Ansatz can be resumed in the following 
theorem. 

Theorem 3.1. The highest weight vector Sl + is an eigenvector of t(X) 

t(X)n + = A {X)Q + (3.1) 

with the corresponding eigenvalue 

A » w - 'w - 2 "'« A » - 1 jh: (t ^7) + 1 m 

a=l \b=jta I a=l 

Furthermore, the action of the B-operators on the highest spin vector tt + yields the Bethe 
vectors 

y M (n)=B M (n)tt+, (3.3) 

such that 

t(\)* u (p) = A M (A; fi)^u(fi), (3.4) 

with the eigenvalues 

A M (A;/x)=pl,(A;/x)-2%l(A;/i) and p M (X; /z) = p(X) - V — ?— , (3.5) 

ox ^ x - m 

i—i 

provided that the Bethe equations are imposed on the parameters \x = {fi±, ...,(1m} 

N M „ 

J2^^-J2^— = 0, i = l,...,M. (3.6) 

a=l & Za j=it & ^ 

Proof. A consequence of $Tty . jSHB|> and (2.20-21) is 

t(x)n+ = (h 2 (x) - 2h'(x))n+ = ( P 2 (x) - 2p'(x))n+ = A {x)n + . 
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The action of i(A) on the Bcthc vcctors ^ m(m) i|3.3|) is evident from lemma l?31 
t(\)* M (ti) = t(\)B M (n)n+ = A (A)tf M (/*) + [t(X),B M (fi)] n+ 

\ ^ MA) ^ 8 \ ^f Af (AU^ ; )) , m . 



A/ 



MIM 

+ ^ 2 E E^-i^'^)*^^) ] Mwm w ), (3.7) 
i=i 

where ^/^(/i.) = B[ k j\fi)Q + and X + (A)f2+ = 0. When the Bethe equations are imposed on 
the parameters /x = {/ii, . . . , /ijf } 

AT M 

^( w;/i W) = £_^__£;^_ =0, i = l,...,M, (3.8) 

o=l Mi Za fri ^ W 

the unwanted terms in Ij3.7|l vanish and hence t(X)^ m (m) = Am(A; fi)^M(^) with 
A M (A; /,) = A (A) - J2 + E a 77T T = ^ ~ ^ ( A; ^ ( 3 ' 9 ) 

□ 

Due to the pole expansion (|2.23|) of t(X) the previous theorem yields the spectrum of the 
Gaudin hamiltoneans. 

Corollary 3.1. The Bethe vectors ^m{p) a-re eigenvectors of the Gaudin hamiltonians 

ffW*MW=4 )$ M(M) (3-10) 

with the eigenvalues 

(3.11) 

Remark 3.1. The fact that the spectrum of the sk-invariant Gaudin model remains is not 
affected by the jordanian twist can also be deduced from the following expression 

t(\)=t(\) Q + 2t(h(\) X+ + 2p^X+(\)) +e(X+) 2 . (3.12) 

Remark 3.2. From the Gaudin realization H2.16(l follows that S,X^ l = lim^oo h(X) and 
together with l|2.3t)|l is straightforward to obtain the action of the global generator X^ on 
the Bethe vectors 

M 
i=l 

Thus the Bethe vectors are annihilated by the global generator X^ 

X+*Af (/*) = (3.14) 
once the Bethe eguations l|3.5|) are imposed. 
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Remark 3.3. Analog ously f rom (|2.16[) follows that — | h g i = liniA^oo X (A) and with l|2.32[l 
and ii. of lemma 12.11 shows that Bethe vectors are not eigenstates of the global generator 

hgl, 

M 

V*m(m) = -2 P [ M) *m(») +2^ M - X) X>M-i(/x (i) ), (3.15) 

i=l 

where p = — h 53^Li and p[ k ^ = (p + k)i = (p + k) ■ ■ ■ (p + l + k — 1) is the Pochhammer 
symbol. 

As it was shown already, the jordanian twist of the s^-invariant Gaudin model preserves 
its spcctrum but the Bcthc vectors are diffcrent, thus their inner products and norms are 
changed also. In order to determinc the inner products and norms of the Bethe vectors it is 
of interest to consider the dual B-operators, obtained by using the dual 

N f X+ f \ 

{x-y{\) = Y J [ T zr7-\K), (3.i6) 



explicitly given by 

M 



B* M {\ U A M ) = H((X-y(\ k ) + (k - 1)0- (3.17) 

fc=i 

As opposed to the s^-invariant case here evidently the dual B-operators do not annihilate 
the highest spin vector il + 

B* M (fi)n+ = z M PM n + . (3.i8) 

Moreover the Bethe vectors are not orthogonal and their norms depend on the twist 
parameter. 

Lemma 3.1. Consider integers M\,M% > 0, M — min{Mi,M2} and complex numbers 
A = {Ai, . . . , Ajv^} and fi = {p\, . . . , (J>m 2 } such that X D fi = 0. Then the inner products 
between Bethe vectors is given by 

M 

n ( k ) J k ) cM 1 +M 2 -2k^ 



fe=i 

Mi M 2 

x Yl Yl (*fc( A ii»---.A i J|*fc(/ij 1 ,... ) A*iJ)o, (3-19) 



where (-|-)o = ("|")|f=Q denotes the corresponding inner product of the sl2-invariant Gaudin 
model 0/. 

Proof. The Bethe vectors are obtained from the action of the B-operators (|2.34[1 on fi + 

M-l M 

*m(/x) = Bm(m)u+ =p M Z M n+ + J2 P ( k M ~ k) e fc E *M- k (» h ,- ■ - 

k=0 h<-<3 M -k 
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where $m(m)o = Bm(p)\^=o^+ are the corresponding Bethe vectors of the s^-invariant 
Gaudin model. Then the inner product of the Bethe vectors of the twisted model are 

Mi-l M 3 -l 
k=0 1=0 

Mi M2 

X £ £ (*M 1 -fc(Ai 1 ,...,Ai Ml _ h )|*M 2 -i(Mji:--- J MjV a -i)}0 

il<---<iMi-fe jl<— <3M 2 -l 

M-l 

.(«) „(") c Mi+M 2 -2n, 



n,, 71», f M i+* _|_ \ " J n > J n > cM 1 +M 2 -2n 
VM X VM 2 K + / j P Mi -nPMi -nS x 

n=l 

Mi M 2 

x £ £ (*»(A 4 i,...,A l J|*„(/j jl ,...,/i 37i )) . 

il<— <in ji< - <j 



The norms of the Bethe vector follow from the previous lemma. 
Corollary 3.2. The norms of Bethe vectors «'m(m) are given by 

M M 
(k) n2 t2(M-k) 



□ 



*M(»)\\ 2 = Y,(p M -k) 2 e iM - k) E <vm^,---,^i*^i,...,/%)>o. (3.20) 



fc=0 ii <•••<«*, 

h<—<jk 



4 Conclusion 

The Gaudin model based on the 5(2 classical r-matrix with jordanian twist is studied. This 
system can be obtained as the semiclassical limit of the XXX spin chain deformed by the 
jordanian twist. Alternatively, applying a certain similarity transformation on the XXZ 
Gaudin model and using the scaling limit, Kulish was able to postulate the spectrum 
and the Bethe states of the system. 

The result of Kulish ^5] that the spectrum of the system coincides with the one of the 
st2-invariant model is demonstrated by full implcmentation of the algebraic Bethe Ansatz. In 
order to construct the Bethe vectors it is necessary to consider the creation operators which 
are not homogeneous polynomials of one of the generators of the corresponding loop algebra. 
However, it was convenient to consider a more general set of operators B^(/j,i, . . . ,/im), 
in order to simplify the calculation of the commutators between the creation operators 
and the generating function of the integrals of motion t(X). These operators are symmetric 
functions of their arguments and they satisfy certain recursive relations. Their commutation 
relations with the generators of the loop algebra are given and they are essential in the 
main step of the algebraic Bethe Ansatz. The creation operators are the particular case 
Bm(mi, ■ ■ ■ j Mm) = B 1 ^ (fj,i, . . . , (j,m)- Thus the commutation relations between the creation 
operators and t(X) are easily calculated. The corresponding Bethe vectors are defined by 
the action of the creation operators on the highest spin vector f2 + and the spectrum of the 
system is determined. However, the Bethe vectors in this case, are not eigenstates of the 
global Cartan generator h g i. 

Some properties of the creation operators are fundamental in calculating the inner prod- 
ucts and the norms of the Bethe states. It is necessary to consider the dual creation op- 
erators -B M (/Ji, . . . , /km) obtained by using the dual Gaudin model. The Bethe vectors 
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*&m(hi , ■ ■ ■ , Mm) are not orthogonal for different M's. Moreover, contrary to the s^-invariant 
case, the generating function of integrals of motion is not hermitian. 

The well known relation ^] between the off-shell Bethe vectors of the Gaudin models 
related to simple Lie algebras and the solutions of Knizhnik-Zamolodchikov equation [231 

d 

k—iI){z 1 , ...,z N ) = H^ip(zu ...,z N ) (4.1) 

OZ a 

where are the Gaudin hamiltoneans (|1.1|) . also holds for the KZ equation related to 

the 5I2 classical r-matrix with the jordanian twist. This relation is obtained by considering 
a Bethe vector ty(jl\z), where the corresponding Bethe equations are not imposed, and the 
integral representation of solutions to the Knizhnik-Zamolodchikov equation, 

tp(z u ...,z N ) = <P(fi\z)y(fiWW (4-2) 

where 4>{fl\z) is a scalar function 

<t>Wz) = ~ Vj) K [[( z a - z b ) « II ||(z a - fi k ) k. (4.3) 

i<j a<b a—l k=l 

The partial derivatives of the scalar factor Ij4.3|l are simply written 

Kd z M=E$(b and Kd w (0) = -2p M ( Mi;i u(%. (4.4) 

It is a simple matter to check that ip given by l|4.2l) satisfies H4.lt . Duc to the Leibniz rulc 
d Za (</)#) = d Za (</>)* + cf>d Za (*) and the residue of (|3~7I) at A = z a 

M 

ffW*(/x)=^ ) *( / x)-2^^( / x)p M (Mi;M W ) (4-5) 

8=1 

where *m (m) = X~(jUi).B$_i(/* W ) fi + and X~(m) = - we have 

M 

k d Za (0*) = (^) + 20 ^ (^)p M M « ) + k a Zo (*) 

1=1 

M 

=ff (o) (^)-«E^ (^M a) )- 



(4.6) 



To obtain the final formula we have used the lemma l2~4l and formulas Q4.4jl . Moreover, 
a closed contour integration of <p^f with respect to fjti, . . . , (j,m will cancel the contribution 
from the sum and therefore ip(zi, . . . , Zn) given by (|4.2|l will satisfy Knizhnik-Zamolodchikov 
equation. 
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A Appendix: proofs of lemmas 

Proof of lemma f2.ll The recurrence relations in ii. and iv. are evident from the definition 

EH 

i. Given a fixed integer k, induction on M is used. Consider M — 2 in definition (|2.24|l : 

B| fe) (/ii,M2) = x-( Pl )x-( P2 ) + (k + i)tx-(m) + k^x-( P2 ) + k(k + i)e, (A.l) 



then using (|2.5I) it is straightforward to check that Bi; (Mi>M2) — B\ J (/t2,/*i). 

Assume B^\ni, . . . , /iat) is symmetric for M — 1 > N. For 1 < i < j < M it is clear 
from ii. that 

^M^(Mi) • • • > Mi, • • • >Mj> • • • > Mm) = #m (Mi, • • • >Mj> • • • > Mi» • • • 1 Mm)- 

The symmetry of B^j\fj,i, . . . ,//m) with respect to /iM-l an d Mm must be shown. To this 
end the recurrence relation in ii. is to be iterated twice and the appropriate terms combined 

Bjj^Oii. • • • > Mm) = SW_ 2 (M1, • • • , MM-2) (X-(^_i) + (M + fc - 2)£) 
x (X-( MM ) + (M + fe - 1)0 = M?-i(Mi, ■ ■ • , Mm-2, Mm) (X-(/i M _i) + (M + fe - 2)0 
+ B^.a^i, ■ • .,HM- 2 )([X-(m-i),X-( m )] + C (X-(/i M -i) + (M + fc - 2)0) 

= #m_i(Mi> • ■ ■ , Mm-2, Mm) (X~(/i M -i) + (M + fc - 2)0 + ^B^^fjn, fi M -2, Mm) 

? (fc) 
'm 



= M?(Mi, • ■ ■ i Mm-2, Mm,Mm-i)- 



iii. Applying the induction on M, for a fbced integer fc. Set M = 2, it is a direct consequence 
of {Ot that 

i?f > ( Ml , /ia) = B^" 1 * ( m , M2 ) + e(B( fe) (Mi) + 0*a)). 
Assume iii. is true for M — 1 then 

^W^B^^'lHE^fM^)) (X-( MM ) + (M + fc-l)0 

M-l 

= 4r^)+^riV M) )+e E (<-i(M w ) + e<- 2 (M ( - M) )) 

4=1 

M-l M 



;.=i 



i=l 



□ 

Proof of lemma 12.21 From definition 12 . II B [ k ^ (u) = X~(fi) + fc£, and commutators (|2.5|l 
it is clear that M = 1 in (2.26-28) is given by 



h(\),B[ k \n) = [h(X),X-0i)] 



B[ k \\)-B\ K '^) 
A — /i 



(fc). 



A — fl 



X-(X),B{ k \n) = [X-(X),X-(n)] =esi*V)-^W(A) 
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The induction method is used to demonstrate lemma Assume that (2.26-28) hold for 
Bj^(fi), M — 1 > N > 1 then, to show that these formulas are valid for M, the recurrence 
relation ii. in lcmma l2~71 is used 

h(\),B$>{ri\ = [MAJ.bSUCm^)] (S 1 (/x M ) + (M + fe-l)0+< ) _ 1 (/x (M) ) 



M-l 



E ^ 1 2 ) (m (i ' m) )[m^)^-(mm)] + <-i(m (m) ) 2 



5i( A) -5i(^m) 
A - m a/ 

i=l i=l ^ 



2 

Mm — M« 



M-l 



Af-l 



■«o« w ) E +q +e E k, m) 

*=i ^ MM v »=1 / 

■f (a fl " (Au f!'- flS ' M + {J»iiti'(M'")3 M ( w m'")) • 



here where used, where appropriate, the induction hypothesis, the commutators [h(X), X (/i)] 
and properties ii. and iv. in lemma l2~71 Thus l|2.26ll is proved. 



X+(\),B$<jj.) = X+(\),B$_ 1 (nM) (B 1 ( m ) + (M + k-l)0 + B^_ 1 (n^) 



3 (fc) 



(k) 



x [X + (X),X~(fiu)] ^ _^ 1 ^+i) (M W ) ^-i(A;M (i - M) ) - 4Af-i(Mi;M (M/) ) 



A - Mi 



M-l r>( fc +l)/r U , , M-l 



A/-1 



/i( A) - %m) 
A — mm 



£X+(A) 



M-l 



- E ^V (l < M) ) 



[/i(A),X-(/x M )]-[/i(Mi),^-(MM)] 



i=l 
M-l 



A - Mi 



E g (m (0 ) ^ M - A; mW - ~ 13 M ^ : M w - - 2 e 



M-l R (fe+1) 



i=l 



Bi^({A}U M (< - J) ) 
(A-Mi)(A-Mj) 



(k+i), (M )jM(\;» iM) )-pM(m;» iM) ) 



A — mm 



+(E<- 1 '(m (i) )^w-<- 1 1 ) (m 
2 ^ ^+;) ({ a}u m (^)) 

^ (A — M») (A ^ Mm) 

7>g-iV V M(A;juW) -fato^ _ 2 f B «-i (W u 

^ ' A /1^ 



i=l 
M 



(A - m)(x - fij) 



+ ^<_ 1 »( M C))I+(A) 



i=l 
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here we used, where appropriate, the commutators [h(X), X (fi)] and [X + (A), X (/*)]. Thus 
(|2~2"7|) is proved. 

'x-(\),B$(n)] = [X-(X),B^_ 1 (^)] (B 1 ( m ) + {M + k-l)0 + B^_ 1 {^) 

x [X-(\), Bl ( m )} = f(M-l)^W_ 1 (^))-^ 1 5W_ 1 ({A}U^))j 

x (B, ( MM ) + (M + fc - 1)0 + £5 M L> (M) )(£i(^Af) - Si(A)) 
M 

=M£<VwE i? M(W u ' u(i) )' 



here we used, where appropriate, the commutators [X (A),X (p)]. 
Proof of lemma 12.41 In the case for M = 1 the proof is straightforward 



□ 



dz a (/i-z a ) 2 djjL\iJLi-z a 2 



Assume (|2.36[) is true for any set /i with |/x| < M complex numbers. Use formula l|2.29[l and 
the induction hypothesis together with 1. of lemma ITTI we can write 



J) 

'd z 



-S M (/x) = (b M -i(m (M) ) (5i(/i M ) + (M - 1)0) 

M-l Q O 

= - E ^7 (^-(MO^-r^)) " 

M-l „ 

=-e|-(^(^X- 1 (^ ) )) 

i— 1 



S M -i(m (M) )^(mm) 
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M 



X-( m )B M ^(n^) + 



[B M ^ M )),X-] 
Mm — 



□ 
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